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We report new data for the topological susceptibility computed on 2+1 flavor dynamical config- 
urations with lattice spacing 0.06 fm, generated with the asqtad action. The topological suscep- 
tibility is computed by HYP smearing and compared with rooted staggered chiral perturbation 
theory as the pion mass goes to zero. At 0.06 fm, the raw data is already quite close to the con- 
tinuum extrapolated values obtained from coarser lattices. These results provide a further test of 
the asqtad action with rooted staggered flavors. 
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1. Introduction 

In 2001 S. Diirr presented [|T]] an analysis of the dependence of the topological susceptibility 
on the pion mass, as measured in then current full QCD simulations. His comparison included 
results using Wilson fermions from CP-PACS, UKQCD, and SESAM/T^L, as well as results using 
thin link staggered fermions from the Pisa group, and by A. Hasenfratz who analyzed MILC and 
Columbia dynamical lattices. The conclusion from these studies was that simulations were not yet 
in agreement with chiral perturbation theory [||] which says that (using f % = 130 MeV) 
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as the pion mass tends to zero. While in most simulations there was a reduction in Z top o as nv\ 
is reduced, contact with the above line was largely absent. This is displayed here in figure 1, 
reproduced from [[I]], with the above expectation (linear in m 2 n ) shown as the left black line, against 
the data. Also shown is the quenched m — > °° expectation as a horizontal line on the right. 
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Figure 1: Comparison of full QCD data for x t0 p ca. 2001, taken from [gj]. 

About the same time, improved actions were seeing a renaissance and have since generated 
striking results in almost all areas of lattice gauge theory. 

In this contribution we present the latest results of the MILC collaboration for the topological 
susceptibility and its mass dependence as m — > 0, using the Asqtad action [Q] at a lattice spacing of 
0.06 fm. With this improved action, HYP smearing, a variance reduction technique to determine 
the topological charge, and careful extrapolation to the continuum limit using rooted staggered 
chiral perturbation theory, we indeed see encouraging agreement with the expectations from QCD. 
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2. Simulations 

By 2003 the MILC Collaboration had lattices at two lattice spacings, a = 0.12 and a = 0.09 fm, 
and a variety of masses with which to investigate the dependence of the topological susceptibility 
on quark (or equivalently pion) mass; these results were presented in [||] and are shown below in 
figure 2. 




Figure 2: MILC data from [f[ (2003). 

This result was quite encouraging, as the expectations of chiral perturbation theory appear to 
be supported by the data. 

In the last 4 years, much larger lattices (up to 48 3 x 144) at a lattice spacing of 0.06 fm have 
been generated and their topological charge analyzed. It is the purpose of this note to update the 
above picture with the new results. In addition to now having three lattice spacings, there have also 
been some theoretical developments in the methods used to extrapolate continuum values which 
we will report as well. 

2.1 Topological charge and susceptibility measurements 

As in we continue to measure the topological charge density F^ V F^ V , using the Boulder 
extended link definition and three HYP smearing sweeps [Sj. 

Typically the topological susceptibility, Z to po = (Q 2 ) /V, is computed by averaging the indi- 
vidual Q 2 from each lattice in the ensemble (V is the lattice volume). For our study, we have used 
a method introduced in [||] which significantly reduces the variance in Q . Since 

(Q 2 ) f 

V = y= J dr(q(r)q(0)) 

we measure the correlator (q(r)q(0)), which we split into a short distance (r < r cut ) part and a long 
distance part (r < r cut ). r cut is ~8-10 in lattice units. 
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At short distance, where the correlator is large, we use the measured points in computing 
(q(r)q(0)}, whereas at long distance since the measured correlator has large variance, we use values 
obtained by a fit to points in the transition region (using the green points in the example shown in 
figure 3). The fit is to a Euclidean scalar propagator, mK\(mr) /4n 2 r, which is the expected long 
distance behavior of (q(r)q(O)). 
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Figure 3: Points used to compute (q(r)q(Q)). Measured points at r < r cut ~ 9a are used. For r > r cut values 
from a fit to those points in green are used. 



2.2 Staggered chiral perturbation theory 

In addition to finer lattices and a variance reduction method for our Z to po measurement, there 
have also been improvements to the expectations from chiral perturbation theory. In particular 
there is now good analytic understanding of the taste splittings induced in operators at finite lattice 
spacing. [^, g|. An example of these taste splittings in the staggered pion multiplet is shown in 
figure 4. The masses of the 16 pions of various tastes are shown as the quark mass is taken to zero. 
At finite lattice spacing the splitting is seen to be roughly constant. As a — > these pions become 
degenerate. The residual t/(l) chiral symmetry of staggered fermions provides one Goldstone pion 
in the multiplet (shown in the figure above with black diamonds), whose mass goes to zero with 
vanishing quark mass, even at nonzero lattice spacing. Typically it has been this pion which is 
studied in lattice QCD, as was the case in our previous topological susceptibility study [||]. 

In however, it was recognized that the staggered pion field which couples to the anomaly, 
and hence should be used for the relevant chiral perturbation theory, is the the taste singlet pion, 
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Figure 4: Taste splittings of the staggered pion taste multiplet as a function of quark mass 



with mass m K j, as opposed to the commonly used Goldstone pion (a pseudoscalar in taste). This 
taste singlet pion is the one displayed with magenta diamonds in figure 4, at the top of the multiplet. 
Following [0, |8|] and starting from the rooted staggered chiral lagrangian 



f 2 11 f 2 m 2 

J? = T -f Tr(d^dU) - ^f- Tr[Jt{U* + U)\ + + £c,^ + . . . 



„2 



where ^0q 7 is an explicit mass term representing the coupling of the anomaly to the taste singlet 
pion field 0o/> Billeter, Detar, and Osborn derive [||] the following dependence of X topo on pion 
masses in the 2+1 flavor case: 



1 + m lj/( 2m ™,i) + 3m l,//( 2m o) 



where m^j is the ss taste singlet pseudoscalar meson mass. This formula interpolates smoothly 
between the m 2 nI —* chiral limit: 

fWnJ 



lim Ztopo 

m->0 



and the quenched limit: 



lim X^-^* 0.06/4 



In the last formula, we use the measured quenched topological susceptibility to set the value of mo, 
and we repeat that m n / here is the mass of the taste-singlet pion. 
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am u j 1 am s 


10/g 2 


lattice spacing 


L 3 x T 


# lats. 


quenched 


8.00 


a = 0.12fm 


20 3 x64 


408 


0.05 / 0.05 


6.85 






425 


0.04 / 0.05 


6.83 






351 


0.03 / 0.05 


6.81 






564 


0.02 / 0.05 


6.79 






484 


0.01 / 0.05 


6.76 






658 


0.007 / 0.05 


6.76 






493 


quenched 


8.40 


a = 0.09 fm 


28 3 x96 


396 


0.031/0.031 


7.18 






496 


0.0124/0.031 


7.11 






527 


0.0062/0.031 


7.09 






592 


0.0072 / 0.018 


7.48 


a = 0.06 fm 


48 3 xl44 


624 


0.0036/0.018 


7.47 






608 



Table 1: Lattices and parameters used in this study 



3. Results 

Putting these developments together, we present our latest results for the 2+1 flavor topolog- 
ical susceptibility. The lattices used for this study have taken more than five years to produce and 
analyze, and are shown in Table 1 . Having Z to po at numerous quark masses and three lattice spac- 
ings, we fit our entire data set to an interpolating function in lattice spacing, a, and taste singlet pion 
mass squared, m 2 nl (the strange taste-singlet mass m s - s ; on these lattices was tuned to be constant): 

y-^(mlj,a) =A + {A l + A 2 a 2 +A 3 a 4 )/"4,/- 

The continuum limit is obtained from the fit by setting a = in this function, and we are left 
with Xtopo {fn 2 nI ) extracted from our data. The result is shown below in figure 5. Measured lattice 
data are shown with blue, red, and magenta symbols, while the continuum limit extrapolation 
function X^o ( m jr /) * s shown with a solid black line. Some representative points along this line are 
shown with error bars reflecting the errors of the continuum extrapolation. Finally, two functions 
representing the chiral perturbation prediction of eq. (2.1) are shown in green: the lower line "L.O. 
2+1+mo" includes the value for mo set by the quenched data, and eq. (2.1) with mo = °° is shown 
labeled "L.O. 2+1" for comparison. 

4. Conclusions 

With the addition of the new a = 0.06 fm data, we see that the topological susceptibility is 
behaving as expected in the m 2 j — > limit of rooted staggered chiral perturbation theory. We find 
it striking that the lightest 0.06 fm datum is almost on the continuum line without extrapolation. 

Finally, we feel that these results lend further credibility to the use of the "fourth root method" 
to simulate single flavors. As mentioned in M. Creutz's talk at this conference, aberrant results 
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Figure 5: Latest MILC results for the topological susceptibility from lattices down to a = 0.06 fm. 

from the fourth root would be expected to arise first in violations of topological quantities and 
correlations, which are quite sensitive to the number of flavors. We see no such violations, and 
indeed only strong support that the simulations are behaving as expected from QCD. 
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